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Many experiments have been done to deter- 
mine the relative strength of different knots, and 
these show that the break in a knotted rope al- 
most invariably occurs at a point just outside the 
'entrance' to the knot 1 . The influence of knots on 
the properties of polymers has become of great 
interest, in part because of their effect on me- 
chanical properties 2 . Knot theory 3,4 applied to 
the topology of macromolecules 5 8 indicates that 
the simple trefoil or 'overhand' knot is likely to 
be present with high probability in any long poly- 
mer strand 9 12 . Fragments of DNA have been ob- 
served to contain such knots in experiments 13,14 
and computer simulations 15 . Here we use ab initio 
computational methods 16 to investigate the effect 
of a trefoil knot on the breaking strength of a 
polymer strand. We find that the knot weakens 
the strand significantly, and that, like a knotted 
rope, it breaks under tension at the entrance to 
the knot. 

Little is known about the structure and properties of 
knots at the atomic level. For example, the minimum 
number of carbon atoms that can be sustained as a tre- 
foil in a polyethylene strand is not well established 5 ' 19 . 
Polyethylene is the simplest polymer and therefore an 
excellent generic system to study the fundamental prop- 
erties of a knotted chain. (Rigorously, a knot is a 
closed loop; but here, as in ref. 1, we apply the term 
to a strand with unlinked ends. As a representative 
polyethylene-like system we chose the linear molecule n- 
decane (CioH^). Empirical models that utilize force con- 
stants for bond stretch, bend and torsion 18-20 , although 
useful for studying bulk properties of chain molecules, 
are not applicable to the case of chain rupture. First- 
principles calculations 16 , on the other hand, have been 
shown to yield satisfactory results for structural and 
mechanical properties of hydrocarbon-based polymeric 
systems 21,22 . The ab initio equilibrium structural pa- 
rameters, after optimization, are in excellent agreement 
with experiment, with errors in bond lengths and angles 
smaller than 1% (ref. 22). We studied the response of the 
n-decane molecule to uniaxial strain by systematically 
increasing the separation L between the terminal carbon 
atoms. For a fixed value of this separation the system 
was heated to room temperature and then cooled via a 
simulated annealing technique. The decane molecule re- 
mained intact for an elongation up to 18.5% beyond its 
equilibrium length. Deviations of the C1-C2 ad C9-C10 



bond lengths from the mean value increase significant 
with the applied strain. Under maximum loading the 
terminal bonds become stretched to ~ 1.8 A compared 
with ~ 1.7A for the others and the terminal C-C-C bond 
angles increase to ~ 135 deE. 
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FIG. 1. Strain energy distribution in a knotted polymer 
strand. Shown are distributions in chains of 35 (a) and 28 
(b) carbon atoms taken from constrained classical MD sim- 
ulations. When the knot is sufficiently tightened, the strain 
energy localizes mostly on the bonds immediately outside its 
entrance points. 

For larger elongation the molecule dissociates into two 
radicals with C10H22 — > CgHig- + CH3-. The bond that 
dissociates involves (randomly) one of the two atoms 
where the tension is applied. We obtain 83 kcal mol -1 
for the dissociation energy for the methyl group, which 
compares favourably with the experimental dissociation 
enthalpy of ~ 87 kcal mol -1 . A similar calculation per- 
formed for n-undecane, C11H24, but with the tension still 
applied to the Ci and C10 atoms, shows that the C9-C10 
bond breaks to yield an ethyl radical. The calculated dis- 
sociation energy in this case is 81 kcal mol -1 compared 
with the experimental enthalpy of ~ 82 kcal mol -1 . As 
in the previous case, the hybridization of the two key 
carbon atoms changes from the initial tetrahedral sp 3 to 
planar sp 2 on formation of the product radicals. 

To generate a starting configuration for the ab ini- 
tio study we first performed classical molecular dynam- 
ics (MD) calculations on a loosely knotted polyethylene- 



like alkane chain, with 7V=144 carbon atoms, in which 
a trefoil was introduced by adding an appropriate set 
of gauche defects 17 . The potential model employed was 
based on the united atom scheme 18,19 but with ab ini- 
tio intra-chain stretch and bend force constants taken 
from the present study. MD calculations were carried 
out on the knotted chain for a sequence of increasing 
(constrained) end-to-end distances, L, which gradually 
resulted in a tighter knot. For any given L, 200 ps of 
room temperature dynamics, followed by 50 ps of slow 
cooling, was utilized to relax the system in the trefoil 
configuration 23 . (This run length seemed appropriate 
given that even the slowest longitudinal phonons can 
propagate along the entire length of the molecule about 
1,000 times during this trajectory.) 
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FIG. 2. Analysis of chain rupture. The breakpoint in 
a 7V=28 trefoil before (a) and after (b) chain dissociation, 
with C22 and C23 and their neighbouring atoms highlighted 
in green and red respectively. Panel c shows the time evo- 
lution of the sum of the respective bond angles, which re- 
veals an asymptotic sp 2 hybridization of the carbons involved 
in the break. The calculations are based on density func- 
tional theory with Becke and Lee- Yang-Parr (BLYP) 24 ' 25 cor- 
rections to the exchange and correlation functionals, respec- 
tively. Electronic spin polarization has been included within 
the local spin-density approximation scheme. A valence elec- 
tron pseudopotential scheme is employed with a plane-wave 
kinetic-energy cutoff at 60 Ry and a T-point sampling of the 
Brillouin zone. The simulation supercell is large enough to 
avoid interactions with periodic images. 

The global (either constrained or unconstrained) min- 
imum corresponds to an unknotted polymer, but the mi- 
gration barrier for the knot to propagate along the chain, 
and thus the disentanglement time, increases dramati- 
cally as the trefoil tightens. Thus even at room tempera- 
ture the system is 'trapped' in the knotted configuration 



for the timescale of the calculations. Under these circum- 
stances, the portions of the chain far from the trefoil, 
although extremely important in disentanglement pro- 
cesses, are likely to have little or no effect on the behav- 
ior of the knot under stress. Accordingly, as the knot 
was tightened, atoms far from it were removed from the 
chain, in order to eventually make the system tractable 
with first-principles methods. 

For each classical MD calculation, we monitored the 
strain distribution along the chain. Even for N = 50, 
the chain with a trefoil stores a rather small amount of 
strain energy. However, when the chain was shortened 
to A^=35 a pronounced inhomogeneity appeared in the 
strain energy distribution (Fig. la). The shape of the 
strain energy distribution is essentially the same, at a 
given L, independent of the details of the simulation, 
such as the initial configuration and temperature or the 
length of the MD simulation or a reasonable change in 
the interaction potentials. (The potentials that we have 
employed in the MD simulations were all fitted to the 
properties of alkane molecules close to their equilibrium 
structures. Thus, in the present case where we are deal- 
ing with highly distorted systems, the model's force con- 
stants are likely to yield only semi-quantitative results.) 
Further shortening of the chain to N=28 produced a tre- 
foil (Fig. lb) with distortions in bond lengths and angles 
still below the critical values that yield bond breaking in 
the linear alkanes. 

The constrained equilibrium positions obtained from 
the 250 ps of classical MD evolution of the C28 n-alkane 
provided the carbon skeleton of the starting configuration 
for the first-principles study of the C28H58 molecule. A 
series of Car-Parrinello MD (CPMD) simulations 16 ' 23 " 25 
were carried out with the separation between Ci and C28 
fixed at distances L between 11.50 and 14.00A In each 
case, the system was allowed to evolve dynamically at 
room temperature, after which it was cooled. Typically, 
during the dynamical evolution individual bonds undergo 
large amplitude thermal fluctuations, which decrease on 
cooling. 

During the room temperature L = 13.50A CPMD run 
for ^=28, dissociation occurred at a bond location, just 
outside the entrance to the knot, well separated from the 
terminal atom where the tension was applied. Figure 2 
shows the region of the knot where the chain actually 
ruptured, both before (Fig. 2a) and after (Fig. 2b) the 
break. The change in hybridization of the C22 and C23 
atoms from sp 3 to sp 2 is clearly evident also in Fig. 2c, 
where we display the time evolution of the system im- 
mediately before and after the break. Figure 3 shows a 
snapshot of the electronic charge density after the break 
occurred. The gap in the charge density between the C22 
and C23 atoms confirms the bond breaking at this po- 
sition. Other simulations carried out with L > 13.50A 
yield breakpoints at one of the two bonds coming out 
from the knot 25 . (Density- functional energy barriers are 



not quantitatively very accurate, but the general qualita- 
tive description of phenomena, as well as the dissociation 
energies involved in such processes, are known to be quite 
reliable.) 




FIG. 3. Schematic electron charge density immediately 
after the break. Carbon and hydrogen atoms are displayed 
with spheres corresponding to their respective covalent radii. 
The green area is a contour plot of the region containing 87% 
of the total electronic charge. A gap in the charge density, and 
thus bond-breaking, is observable between the two highlighted 
C atoms. 

The first-principles CPMD calculations provide the to- 
tal strain energy but not its profile along the chain. The 
latter can be obtained by making use of an ab initio force 
constants model, which is able to reproduce the total 
strain energy to ~ 20%. Figure 4 shows the evolution 
of the strain energy distribution during structural relax- 
ation of a iV=30 trefoil just before breaking occurs. The 
strain energy is mainly localized in the two symmetric 
bonds that are outside the entrance to the knot. Our 
calculations suggest that 23 carbon atoms form the tight- 
est knot that can be sustained in a polyethylene strand 
without it breaking. 

The strain energy stored in the knot at breaking point 
is 12.7 kcal mol -1 per C-C bond, which is considerably 
smaller than the value 16.2 kcal mol -1 for the linear un- 
knotted case. Thus, we find that the presence of the knot 
has significantly weakened the strand in which it is tied. 

While the dynamical evolution of the present con- 
strained ab initio MD simulations was sufficient to ob- 
serve bond breaking, no attempt was made to allow for 
recombination of the resulting radicals or reactions of the 
radicals with other parts of the chain. The study of these 
effects as well as the role of chain branching and the in- 
fluence of neighbouring chains is left for future research. 
These factors will likely provide a deeper understanding 
of how the interplay between inter- and intra-molecular 
effects contributes to the mechanical properties of real 
polymer samples 2 . 



Note added in proof: Arai et al. 26 have recently 
reported the knotting of actin filaments and DNA 
molecules using optical tweezers. They find that the 
breaking stress for actin is significantly lower than that 
of the unknotted filaments, and that the breakage point 
is at the entrance to the knot, as our calculations predict. 




Bond number 

FIG. 4. Relaxation of the strain energy distribution. 
Shown is the distribution along a C30HB2 chain containing 
a trefoil during CPMD energy minimization at T = OK. The 
t = configuration is an average of the equilibrium positions 
obtained along classical MD simulations performed with dif- 
ferent initial states and temperatures. The segments outside 
the knot and its central portion relax, while stress tends to 
concentrate mostly on the two entrance bonds (see Fig.l). 
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